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First detection of a low-mass stellar halo around the young open 
cluster 77 Chamaeleontis 
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ABSTRACT 

We have identified several Hthium-rich low-mass (0.08 < M < 0.3 Mq) stars within 5.5 deg 
of the young open cluster r] Chamaeleontis, nearly four times the radius of previous search 
efforts. Of these stars we propose 4 new probable cluster members, and 3 possible mem- 
bers requiring further investigation. These findings are consistent with a dynamical origin for 
the current configuration of the cluster, without the need to invoke an abnormal Initial Mass 
Function deficient in low-mass objects. Candidates were selected on the basis of DENIS and 
2MASS photometry, NOMAD astrometry and extensive follow-up spectroscopy. 
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1 INTRODUCTION 

, The recently-discovered open cluster 77 Chamaeleontis is one of the 
' closest {d ^ 94 pc) and youngest (t ^ 8 Myr) stellar associa- 
tions in the Solar neighbourhood. A census of its stellar population 
currently stands at 18 systems, covering spectral types B8-M5.5. 
At high and intermediate masses the cluster Initial Mass Function 
(IMF) follows that of other star-forming regions and young stellar 
groups, but there is a clear deficit of members at masses <0. 15 Mq . 
Comparing the obse rved mass function of the cluster to the Orion 
Trapezeium Cluster JLvo et al.l (l2004h predict that an additional 20 
stars and brown dwarfs in the mass range 0.025 < M < 0.15 
Mq remain to be discovered. Efforts to observe this hitherto un- 
seen population have failed to find any additional members at ei- 
ther larger radii from the cluster core lLuhmanll2004 to 1.5 deg, 
4 times the r adius of known membership) or to low masses in the 
cluster core (ILvo et al.ll200d to ^-^13 Mjup). Failure to find these 
low-mass members raises a fundamental question: has the cluster's 
evolution been driven by dynamical interactions which dispersed 
the stars into a diffuse halo at even larger radii, or does 77 Cha 
possess an abnormally top-heavy IMF deficient in low-mass ob- 
jects? The latter result would seemingly be at odds with the grow- 
ing body of evidence that suggests the IMF is universal and inde- 
pendent of initial star-forming conditions (for a detailed review see 
iBastian. Covev & MeveJl201Qh. 



Moraux. Lawson & Clarkd (l2007l) have attempted to model 



the observed properties of 77 Cha using A^-body simulations of the 
cluster's dynamical evolution starting with standard initial condi- 
tions. They are able to replicate the current configuration of the 



cluster assuming a log-normal IMF and 30-70 initial members. 
This suggests the deficit of low-mass objects seen in the present 
day cluster may not be due to a peculiar IMF but to dynamical evo- 
lution. Their simulations predict there should exist a diffuse halo 
of cluster ejectees beyond the radius currently surveyed. Such co- 
moving objects will have similar proper motions and distances as 
the cluster proper. In this letter we present the results of our search 
for this putative halo of low-mass stars surrounding 77 Cha. 



2 CANDIDATE SELECTION 

2.1 Photometry 

Using the ToPCATlj Virtual Observatory tool we crossmatched DE- 
NIS and 2MASS survey photometry within a 5.5 deg radius sur- 
rounding 77 Cha. Candidates were selected by forming the zdenis 
versus (i — J2mass) colour-magniude diagram for the 1.2 x 10^ 
sources in the field. As seen in Figure [T] the empirical cluster 
isochrone is well defined in this diagram through the K and M spec- 
tral types. Due to its youth, proximity to Earth and lack of any sub- 
stantial reddening, the cluster isochrone sits well above the vast ma- 
jority of stars in the field. To account for unresolved binarity, red- 
dening, distance variations and any non-linearities in the isochrone 
we initially selected photometric candidates lying ±1.5 mag from 
a linear fit to all known late- type members. Because we are inter- 
ested in the low-mass population we restricted our candidates to 
{i — J) > 1.5, corresponding to spectral types ^M3 and later, 
and (J — ii^)2MASS < 0.8, to minimize contamination from back- 
ground giants. Excluding stars with poor 2MASS photometry (14), 
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Figure 1. DENIS/2MASS colour-magnitude diagram for a 5.5 deg ra- 
dius around 77 Cha. Contours show the cumulative totals of stars enclosed. 
We select candidates (filled circles) ±1.5 mag from the empirical cluster 
isochrone (solid line) also having (i — J) > 1.5. Known KM-type mem- 
bers are shown as filled stars. Proper motion candidates are denoted by open 
circles. The intermediate-gravity stars described in ^3.3l are labelled. 
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tain proper motion measurements explicitly set to zero (D. Monet, 
private communication). All 18 77 Cha members were found in NO- 
MAD, 14 of which have non-zero proper motions. The four out- 
standing members are located close to bright stars or diffraction 
spikes, making proper motion measurements difficult. 

For a given UVW space motion the resultant proper mo- 
tion vector (and radial velocity) depends on both sky position 
and distance. Given the large angular extent of the survey area, 
this can have a substantial effect on the expected proper motions 
and radial velocities of our candidates. Using the canonical clus- 
ter space motion, (U V, W) = (-11.8, -19.1, -10.5) km s"^ 
(iMamaiek et al.ll200Qi) we calculated the expected proper motion 
of each candidate. From the uncertainties in the NOMAD proper 
motions we then calculated the number of sigma from the expected 
proper motion in each component. We selected as proper motion 
candidates the 14 stars that lay within ±3a of the mean of the clus- 
ter stars. All 14 known members fell within the 3 a selection box. 



3 FOLLOW-UP SPECTROSCOPY 

3.1 Low-resolution spectroscopy 

We obtained low-resolution spectra of our proper motion candi- 
dates and a selection of photometric candidates using the Double 
Beam Spectrograph (DBS) on the ANU 2.3 m telescope at Siding 
Spring Observatory during 2009 March. With the 316 lines mm~^ 
red grating our spectra cover the wavelength range 5920-9750 A at 
a resolution of 5 A (1.9 A pixel" ^). On each night we also observed 
flux standards, smooth spectrum stars for telluric co r rection and a 
selection of Gliese dwarf standards from ILvo et al.l (|2004 l2008h 
for spectral typing and gravity estimation. Of the 14 proper motion 
candidates, 1 1 were observed in this mode. The remaining three are 
faint and two have DENIS i magnitudes that fall more than 1 mag 
below the cluster isochrone and are not considered likely mem- 
bers. We also observed 18 additional photometric candidates that 
lay close to the cluster isochrone and had null proper motions, or a 
proper motion just outside the 3 a selection box. 



Figure 2. Distribution of candidates on the sky. Symbols a re as for F i gured 
The area enclosed by the dotted line is that searched bv iLuhmanI (l2004l) . 
the dashed line is the 5.5 deg radius of this study. Contours correspond to 
10 MJy sr-i IRAS 100 M-m flux. 



known members (5) and previously confirmed field stars (3) left 
81 photometric candidates. Figure |2] shows the distribution of the 
candidates on the sky. Just one candidate l ies within 1.5 deg of the 
cluster, the radius to which oihmani (|2004|) searched using a differ- 
ent method and reported no new members. 



2.2 Proper Motions 

Kinematics have been vital in confirming membership of the nu- 
merous young stellar associations in the Solar neighborhood. In 
the absence of radial velocities and accurate distances however, 
proper motions are the sole means of investigating the kinematics 
of our stars. We therefore cross matched the 81 candidates against 
the Naval Observatory Merged Astrometric Catalogue (NOMAD). 
Only one of the photometric candidates was not found in NOMAD, 
this is the second faintest star in the sample. Eighteen candidates 
returned a null proper motion. These correspond to highly uncer- 



3.2 Spectral types 

To garner spectral types we first performed synthetic photometry 
on the flux-calibrated DBS spectra to calculate broadband (R — I) 
colours. Unfortunately, the Gliese standards we intended to use to 
calibrate the colour transformations are much brighter than the can- 
didates and suffer from severe fringing at wavelengths greater than 
^8000 A, probably as a result of scattered light inside the spec- 
trograph. Using the Gliese dwarfs to calibrate the candidate pho- 
tometry results in a systematic 0.05-0.15 mag over-estimation of 
(R — I) for the known 77 Cha members. We therefore calibrated 
the photometry of the Gliese dwarfs and candidates separately. The 
instrumental (R — I) colours of the two groups were transformed 
onto the Cousins system using published photometry for the Gliese 
dwarfs ([Bes sell 1990) and the 14 late-type members of 77 Cha from 
ILvo et al.! ([2004, 2008), respectively. In both cases the transfor- 
mation was well-fitted by a low-order polynomial. The ILvo et alj 
r] Cha colours are calibrated using the same Gliese dwarfs and 
Bessell (1990) colours we used, so the two calibrations are equiva- 
lent. 

Spectral types of the candidates were finally obtained by con- 
verting the (R — I) colou r to a s p ectral type using the M -dwarf 
transformation of iBesselll (Il99lh . ILvo et all (|2004 l2008h have 
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shown that the (i? — /) colours (and various narrowband spectral 
indices) of intermediate-age pre-main sequence (PMS) stars are 
indistinguishable from Gyr-old field dwarfs and that they closely 
follow a dwarf temperature (spectral type) sequence. The (R — I) 
derived spectral types of the Gliese M-dwarfs and known 77 Cha 
members agree with published values at the 0.1-0.2 subtype level. 
No attempt has been made to correct the photometry for the effects 
of reddening - we assume it to be minimal at these distances. The 
maximum reddening from dust maps is E(B — V) = 0.5 mag, al- 
though the majority of this material is probably background to the 
candidates. We estimate 0.1 mag as a more likely upper limit. To 
include the effects of unknown reddening we adopt 0.5 subtypes as 
the uncertainty on our spectral type measurements. 

3.3 Gravity estimation and Ha emission 

Pre-main sequence stars are in the process of contraction to their 
eventual main sequence radii. As a result they have slightly lower 
surface gravities than their main sequence counterparts (but still 
much closer to dwarfs than giants). The neutral alkali metal lines in 
our spectra, such as the K I (7665/7699 A) and Na I (8183/8195 A) 
doublets, are highly sensitive to gravity. Later than mid-M spec- 
tral types there is a marked difference in Na I absorption between 
M-dwarfs, PMS stars and M-giants. Many authors have used the 
strength of the Na I doublet to gauge the surface graviti es of late- 
type PMS candidates and assign coarse age rankings (e.g. Lyo et al. 
I2OO4 . 2008: Lawson et al. 2009) . We adopt the spectral index de- 
fined by ISlesnick et al.l (l2006h to measure the strength of Na I ab- 
sorption. It is formed by integrating the spectrum in two 30 A 
bands, one on the doublet, the other sampling the adjacent pseudo- 
continuum. Figure [3] shows the Na I index versus (R — I) colour 
for our DBS targets, known 77 Cha members, the Gliese field dwarfs 
and KM giants from the Pickles library. At spectral types later than 
M3 the locus of 77 Cha members clearly diverges from the dwarf 
sequence and assumes an intermediate gravity between dwarfs and 
giants. Several candidates occupy the same region of the diagram - 
these 6 stars are likely PMS stars of ages 5-10 Myr. 

All of the intermediate-gravity stars exhibit Ha emission, as 
do all the late-type members of 77 Cha. This is a necessary, but in- 
sufficient indicator of youth because of the potential contamination 
from older, flaring M-dwarfs. When combined with gravity infor- 
mation however it is a powerful discriminant. The six intermediate- 
gravity stars and their particulars are presented in Table [T] 
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Figure 3. Gravity-temperature diagram for our DBS targets. Thin lines are 
fits to the observed K and M-type dwarfs (squares) and Pickles giants (tri- 
angles). The thick line is a fit to the late-type cluster members (filled stars). 
Candidates with Ho; emission are given by filled circles, those without are 
crosses. Proper motion candidates are denoted by open circles. 



a heliocentric frame). Stars were observed on multiple nights and 
we took the velocity as the mean of these values. Multiple measure- 
ments and cross-correlations between standards show we are easily 
able to attain a radial velocity precision of 2 km s~^ from the in- 
strument. Li I and Ha equivalent widths and the radial velocities 
derived from the WiFeS spectra are given in Table [T] 



4 CLUSTER MEMBERSHIP 

4.1 Lithium equivalent widths 

Lithium is very easily destroyed in stellar atmospheres, hence the 
detection of lithium absorption in a low-mass star is a strong in- 
dicator of youth. The equivalent width of the Li I 6708 A feature 
is highly dependent on both age and temperature. iMentuch et alJ 
(2008) present the distribution of Li I EWs for a variety of nearby 
stellar associations. For the spectral types under consideration here 
all but one of the measured Li I EWs in Table [T] are consistent with 
an 8-10 Myr population like 77 Cha. Based on its smaller Li I equiv- 
alent width and poor CMD placement we are able to immediately 
rule out 2MASS 0942-7727 as an ry Cha member. 



3.4 Medium-resolution spectroscopy 

To derive radial velocities and check for the presence of lithium ab- 
sorption at 6708 A, we obtained medium-resolution (R ^ 7000) 
spectra of the stars in Table \T\ with the WiFeS instrument at the 
ANU 2.3 m during 2010 January- April. We also observed many of 
the other candidates with congruent pho tometry, proper motions or 
Ha emission. WiFeS (DopitaetaLzOOT) is a double-beam, image- 
slicing integral field spectrograph, and provides a 25x38 arcsec 
field with 0.5^^ pixels. The red i?7000 grating provides wavelength 
coverage from 5300-7060 A, at a resolution of 1 A (0.44 A px~^). 
For deriving radial velocities, we took NeAr arc frames immedi- 
ately after each exposure and also observed 5-7 mid-M radial ve- 
locity standards each night. Only the 6 intermediate-gravity stars 
in Table [1] were found to have lithium in their spectra. To obtain a 
radial velocity each spectrum was cross-correlated against all the 
standards from that night over the region 6000-6500 A, and the ve- 
locity taken as the mean of those measurements (after correcting to 



4.2 Other young stars in the region 

The region s urrounding 77 Cha is rich in pre-main sequence stars of 
various ages. lCovino et alj (119971) present radial velocities, spectral 
types and lithium equivalent widths from high-resolution spectra of 
ROS AT detections in the region. Four of their stars lie in the eastern 
half of our survey area and have spectral types and lithium equiv- 
alent widths similar to the late-type population of 77 Cha. These 
are shown in Table[T]with measurements derived from WiFeS spec- 
tra. Spectral types are from Covino et al. (1997). Such an X-ray 
selected sample is probably incomplete, as the ROSAT flux limit 
for a ^10 Myr PMS star at the cluster distance is around M2. 



4.3 Dynamics 

If the remaining candidates are in fact ejected members of 77 Cha 
we do not expect their space motions to be identical to that of the 
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Table 1. Candidate 77 Chamaeleontis members in the surveyed region 



2MASS 
identifier 



^DENIS 



J) Sp. Type. 
±0.5 



HaEW 

[±lA] 



LilEW 
[± 0.05 mA] 



Vr 
[±2kms-i] 



Ma M5 

[± 9 mas yr~^] 



Membership? 
(reason) 



J08014860-8058052 


13.17 


1.88 


M4.4 


-20<EW<-7 


0.55 


19.3 


-28 


+39 


possible (binary?) 


J08202975-8003259 


12.96 


1.65 


M4.3 


-40<EW<-15 


0.60 


18.4 


-33 


+30 


probable 


J09053087-8 134572 


14.16 


2.01 


M4.9 


-9 


0.65 


16.2 


-35 


+27 


probable 


J09 133435-7550099 


13.82 


1.97 


M4.8 


-8 


0.70 


12.9 


-40 


+34 


probable 


J09424157-7727130 


14.11 


1.73 


M4.6 


-7 


0.35 


18.3 


-16 


+23 


no (Li, CMD) 


J09553919-7622119 


13.12 


1.85 


M4.1 


-5 


0.60 


8.8 


-30 


+ 19 


possible (dynamics) 


Covinoetal(1997)RC 


^SAT candidates: 
11.55 1.45 


M3 


-2 


0.55 


18.7 


-35 


+31 




RX J0902.9-7759 


probable 


RXJ0915.5-7608 


10.44 


1.11 


K7 


-1 


0.55 


22.0 


-32 


+ 16 


no (too far) 


RX J0942.7-7726 


11.49 


1.12 


MO 


-3 


0.45 


19.1 


-24 


+ 17 


no (Li, CMD) 


RXJ1005.3-7749 


11.17 


1.36 


Ml 


-4 


0.60 


17.7 


-30 


+ 17 


possible (dynamics) 
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Figure 4. The results of dynamical simulations for the 5 Li-rich members and 4 PMS stars from lCovino et al.l (Il997h . The colour sc ale indicates the diffe rence 
between the observed and expected space motion at that distance. The dashed line denotes the revised 94.3 pc cluster distance from lvan LeeuwenI l2007h . 



cluster proper. They will have each been imparted some ejection 
velocity which acts over time to disperse the star away from the 
cluster core. To determine the possible epoch and magnitude of this 
impulse we modeled the space motion of the candidates as a func- 
tion of ejection time and current distance, between — 10 and Myr, 
and 50 < d < 150 pc. Figured shows the results of such simula- 
tions. At each ejection time the cluster was first linearly backtraced 
to its heliocentric position at that epoch. For each distance we then 
calculate the current Galactic location of the candidate and the cor- 
responding UVW required to move it there since ejection. This 
is compared to the observed UVW derived from proper motion 
and radial velocity measurements. The colour-map in FigurelHindi- 
cates the vector magnitude of the difference between the two space 
motions. Empty regions in the diagrams have velocity differences 
larger than those expected from the error in space motion at that dis- 
tance. A similar map can be made for ejection speed. In these sim- 
ulations we use the new cluster distance of 94.3 pc a nd proper mo- 
tion vector from the revised Hipparcos astrometry of Ivan LeeuwenI 
(12007) and the improved radial velocity, Vr = 18.3 ± 0.1 km s"^ 
(A. Brandeker, private communication). Proper motions from the 



PPMXL catalogue fe oeser et al.ll2010h are also now available. We 
use these here and in Table [T]in preference to the original NOMAD 
values. Although they quote larger errors than NOMAD, PPMXL 
is better able to replicate the mean proper motion of known 77 Cha 
members. As a check we have repeated the proper motion anaysis 
of ^2.21 using the PPMXL values and find no change in the final 
candidate list. 



The kinematic distances estimated from Figure |4] can be 
checked against those expected from a star's position in the clus- 
ter CMD. Such a diagram is shown in Figure [S] where we separate 
the cluster members into single (or high mass-ratio binary) stars 
and binary systems. Equal mass binaries should lie 0.75 mag above 
the single star locus, corresponding to the dashed line. There is ex- 
cellent overall agreement between the kinematic and photometric 
distances of the candidates. From comparison of Figures |4] and [5] 
we propose 4 probable and 3 possible new 77 Cha members. 



First detection of a low -mass stellar halo around rj Cha 




1.0 1.2 1.4 1.6 1.8 2.0 

^DENIS~'^2MASS 



2.2 



Figure 5. The late-type 77 Cha cluster sequence, for single (gray stars, solid 
fitted line) and known binaries (filled stars, dashed line). The arrow indi- 
cates a reddening of E(B — V) = 0.1 mag in this colour-magnitude plane. 
The vertical line is the blue limit of our photometric/proper motion survey. 



RX J 1005. 3. Figure [5] suggests a distance slightly in front of 
the cluster, but the best fitting kinematic distance is 95 pc, corre- 
sponding to an ejection of 4-5 km s~^ 1.5 Myr ago. It lies at the 
edge of the 100 /xm IRAS flux in Figure |2] so it is possible that the 
star may be reddened slightly (see the arrow in Figure |5]). 



4.6 Non-members 

RX J0915.5. The dynamical simulation for this star suggests a dis- 
tance of 1 10 pc, whereas its position in Figure [^indicates it should 
lie at a similar distance to the cluster, some 15 pc closer. 

RX J0942.7. Not a likely 77 Cha member due to its position a 
magnitude below the cluster isochrone and lower lithium equivalent 
width. Even at the 120 pc given by the simulations it still lies 0.5 
mag below the single-star isochrone. It may be a member of the 
older, more distant Lower-Centaurus-Crux OB association. 



4.4 Probable members 

2MASS 0820-8003. This candidate is a highly probable member 
based on its observed space motion. It also has the closest angular 
separation from the cluster core of any candidate (1.44 deg). The 
results of the simulations show that it need only be ejected from the 
cluster at 1-2 km s~ ^ in the last 4 Myr ago to move to its present lo- 
cation. Multiple WiFeS observations show the star possesses strong 
and highly variable Ha emission. 

2MASS 0905-8134. this star shows excellent agreement be- 
tween its kinematic distance of ^-^86 pc and slightly elevated posi- 
tion above the 77 Cha CMD. The simulations suggest the star was 
ejected from the cluster 6.25 Myr ago at a speed of 1.5 km s~^. 

2MASS 0913-7550. also shows excellent agreement between 
kinematic and photometric distances. The star was ejected 3.5 Myr 
ago at a speed of 6 km s~^ and currently lies in front of the cluster 
at a distance of 73 pc. Stellar encounters would be so rare and weak 
at this epoch that the only possible mechanism for such dynamic is 
from the break-up of an unstable multiple system, perhaps one of 
the currently seen binaries in the cluster 

RX J0902.9. This M3 star from ICovino et all (Il997h lies just 
outside the 1.5 deg surveyed bv lLuhmanI (|2004|) . It has a kinematic 
distance of ^85 pc and would have been ejected at 1-2 km s~^ 
early in the cluster's dynamical evolution. 



4.5 Possible members 

2MASS 0801-8058. This star is a probable close binary, based on 
its location in Figure [5] and 8 velocity measurements during 2010 
January- April. The velocity given in Table [T] is the mean of these 
measurements, and should be close to the systemic velocity. The 
dynamical simulations suggest it was ejected from the cluster at a 
speed of <2 km s~^ early in the cluster's evolution. Even taking 
into account the apparent binarity, its position in Figure [5] is more 
consistent with a distance at or slightly behind the cluster, rather 
than the 85 pc given by Figure |5] 

2MASS 0955-7622. Also a possible binary, from its elevated 
position in Figure \5\ and low velocity. However, ten RV measure- 
ments between 2010 January- April do not show any significant 
variation outside instrument errors. If the simulations are correct, 
it would have been ejected from the cluster at over 10 km s~^ very 
late in the cluster's history. Even so, its current distance would still 
not be close enough to match its position in Figure [5] 



5 CONCLUSIONS 

We have discovered several new probable members of 77 Cha up to 
5 deg from the cluster core and several possible members await- 
ing further confirmation. Are these findings consistent with a dy- 
namical origin for the paucity of low-mass objects observed in the 
cluster? 

To ans wer this we consider the radial distribution of ejectees 
presented in lMoraux et alj (12007 ). Integrating the distribution out to 
9 pc (5.5 deg) we could expect to find up to 6 stars across all masses 
in the survey area. Our blue photometry limit of i — J > 1.5 cor- 
responds to an approximate spectral type of M3, whereas our spec- 
troscopic campaign is complete to i — J < 2, approximately M5. 
Transforming these spectral types to temperatures usi ng the relation 
oflBesselll (Il991i) and then to masses from an 8 Mvr lBaraffe et alJ 
(Il998h isochrone gives a n approximate surveyed mass range of 
0.08 < M < 0.3 M0. ' Moraux et all find the mass distribution 
of ejectees is ro ughly constant with radius and consistent with the 
IChabrieJ (l2003h input IMF. We therefore integrate this IMF over 
the above mass range to find the fraction of stars at those masses, 
this is ^40%. Hence we can expect to find 2-3 bona fide 77 Cha 
members within the surveyed area and mass range. Our discov- 
ery of 3 probable 2MASS/DENIS members is consistent with this 
prediction. We can therefore conclude that dynamical evolution is 
solely responsible for the current configuration of 77 Cha and it is 
not necessary to invoke an IMF deficient in low-mass objects. 
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